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Abstract 


Let M be a compact Riemannian manifold with smooth boundary, 
and let /?(A) be the Dirichlet-to-Neumann operator at frequency A. 
The semiclassical Dirichlet-to-Neumann operator R sc i(A) is defined 
to be A~ 1 /?(A). We obtain a leading asymptotic for the spectral 
counting function for R SC \(X) in an interval [a i, 32 ) as A —> 00 , under 
the assumption that the measure of periodic billiards on T* M is zero. 
The asymptotic takes the form 

N(A; ai, 32 ) = (^( 32 ) — k(3i)) vo\'(dM)\ d ~ 1 + o(A 0,_1 ), 
where n(a) is given explicitly by 


«(a) 


Ud-l 
( 2n) d ~ 1 



di1 


-1 + H(a)( l + a 2 )^'^. 
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1 Introduction 

Let M be a Riemannian manifold with boundary. The Dirichlet-to-Neumann 
operator is a family of operators defined on L 2 (dM) depending on the 
parameter A > 0. It is defined as follows: given f E L 2 (dM), we solve the 
equation (if possible) 

(1.1) (A — A 2 )u = 0 on M, u\ 9M — f. 

Then the Dirichlet-to-Neumann operator at frequency A is the map 

(1.2) F(A) : f ^ — 1 9 u u\qm- 

Here d v is the interior unit normal derivative, and A is the positive Laplacian 
on M. 

It is well known that /?(A) is a self-adjoint, semi-bounded from below 
pseudodifferential operator of order 1 on L 2 (dM ), with domain /-/ 1 (9M). It 
therefore has discrete spectrum accumulating only at +oo. The Dirichlet- 
to-Neumann operator and closely related operators are important in a 
number of areas of mathematical analysis including inverse problems (such 
as Calderon’s problem [Cal]), domain decomposition problems (such as the 
determinant gluing formula of Burghelea-Friedlander-Kappeler [BFK]), and 
spectral asymptotics (see e.g. [Fried]). 

In this paper, we are interested in the spectral asymptotics of R{ A) in the 
high-frequency limit, A —> oo. Let us recall standard spectral asymptotics 
for elliptic differential operators, for simplicity in the simplest case of a 
positive self-adjoint second order scalar operator. Suppose that Q is such an 
operator on a manifold M of dimension d, with principal symbol q. Then in 
the case that M is closed, we have an asymptotic for the number N(A) of 
eigenvalues of Q (counted with multiplicity) less than A 2 : 

(1.3) N(A) = (27r)” d vol{(x,0 E T*M \ q(x,£) < A 2 } + 0(A d “ 1 ) 

= (^) d v°l{(*,0 € T-M I ,(*.£) < 1} + OlA^ 1 ), 

where the equality of the two expressions on the RHS is a simple consequence 
of the homogeneity of q. Moreover, if the set of periodic geodesics has 
measure zero, then there is a two-term expansion of the form 

(^) dvo| {(*.0 e T * M I <7(*.0 < 1} + I sub(Q) + o(A d ” 1 ) 

2vr (27r) d J {q=1} 
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where sub(Q) is the subprincipal symbol of Q [DuGu]. This was generalised 
to the case of manifolds with boundary by the second author [Ivr80]. For 
simplicity we state the result in the case that Q = A is the (positive) metric 
Laplacian, which satisfies sub(A) = 0. Then A is a self-adjoint operator 
under either Dirichlet (—) or Neumann (+) boundary conditions, and if 
the set of periodic generalised bicharacteristics has measure zero, we get a 
two-term expansion for Na(A) of the form 

(1.4) ( — ) d voie*/W± -(^-) d ~ 1 vo\B*dM + o(\ d ~ 1 ). 

' 2.71' 4 x 2t[' 

These statements can be generalized to the semiclassical setting. Consider 
a classical Schrodinger operator on M, P = h 2 A + l/(x) — 1, where h > 0 
is a small parameter (“Planck’s constant”) and V is a smooth real-valued 
function. We consider the asymptotic behaviour N^(P) of the number of 
negative eigenvalues of P as h —y 0. This is equivalent to the problem above 
if h = A -1 and 1/ is identically zero. Define p(x, £) to be the semiclassical 
symbol of P, i.e. p — + V'(x) — 1- Then, if M is closed, under the 

assumption that the measure of periodic bicharacteristics of P is zero in 
T*M, and that 0 is a regular value for p, we have 

(1.5) N h (P) = (2vr/))~ cy vol{(x,0 e T*M \ p(x,£) < 0} + 0(/) 1 “ d ). 

Moreover, for manifolds with boundary, we have an analogue of (1.4): under 
either Dirichlet (—) or Neumann (+) boundary conditions, if the set of 
periodic generalised bicharacteristics has measure zero, we get a two-term 
expansion for N h (P) (where here we understand the self-adjoint realization 
of P with either Dirichlet or Neumann boundary condition) of the form 

(1.6) (2vr/?)" d vol{(x,0 e T*M \ p(x, 0 < 0} ± ^vrh) 1 ^ vol PL + o(/7 1 ~ rf ), 

where H C T*(dM) is the hyperbolic region in the boundary, that is, the 
projection of the set {(x,£) | p(x, 0 < 0} D Tq M M to T*0M. 

From the semiclassical point of view, since P(A) is a first order operator, 
it makes sense to consider P sc i(A) A _1 P(A) (for A > 0), which we call the 
semiclassical Dirichlet-to-Neumann operator. Like P(A), it is a self-adjoint, 
semi-bounded from below operator on L 2 (0M), with discrete spectrum 
accumulating only at +oo. The goal of this paper is to investigate the 
spectral asymptotics of P sc i(A), that is, the asymptotics of 

(1.7) N( A; ai, a 2 ) := #{/3 : f3 is an eigenvalue of P sc i(A), 3i < fd < a 2 }, 
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the number of eigenvalues of R sc i(A) in the interval [ai, a 2 ), as A —> oo. 

Both R( A) and /? sc i(A) have the disadvantage that they are undefined 
whenever A 2 is a Dirichief eigenvalue, since then (1.1) is not solvable for 
arbitrary f G H 1 (M). Indeed, when A 2 is a Dirichlet eigenvalue, a necessary 
condition for solvability of (1.1) is that f is orthogonal to the normal 
derivatives of Dirichlet eigenfunctions at frequency A. To overcome this 
issue, we introduce the Cayley transform of R sc |(A): we define 

(1.8) C(A) = (/? SC ,(A) - /)(/?sd(A) + O' 1 - 

This family of operators is related to impedance boundary conditions: we 
have C(X)f — g if and only if there is a function u on M satisfying 

(1.9) (A — X 2 )u — 0 
and 

(1.10) ! 2 ^(X~ 1 d u u - iu) — f, ^(A _1 <9 U u + iu) = g. 

As observed in [BH], C(A) is a well-defined analytic family of operators for 
A in a neighbourhood of the positive real axis, which is unitary on the real 
axis. In particular, it is well-defined even when A 2 is a Dirichlet eigenvalue 
of the Laplacian on M. As a unitary operator, C(A), A > 0, has its spectrum 
on the unit circle, and as R sc |(A) has discrete spectrum accumulating only at 
oo, it follows that the spectrum of C(A) is discrete on the unit circle except 
at the point 1. Our question can be formulated in terms of C(A): given two 
angles 6\, 9 2 satisfying 0 < 9\ < d 2 < 2i r, what is the leading asymptotic for 

(1.11) N(A;0i,0 2 ) 4f{ e ' d '■ e '° is an eigenvalue of C(A), 9 1 < 9 < 9 2 } 

the number of eigenvalues of C(A) in the interval {e' e : 9 e \9\, 9 2 )} of the 
unit circle, as A —» oo. Clearly, we have 

(1.12) N(A; 9i, 9 2 ) = N(A; a lT a 2 ), where e ldj = — —i.e. 3j — — cot(^). 

To answer this question we relate it to a standard semiclassical eigenvalue 
counting problem on M. To state the next result, we first define the self- 
adjoint operator P 3 h on L 2 (M) by 

(1.13) D(P 3i/i ) = {u G H 2 (M) : ( hd v + a)u = 0 at dM}, 

(1.14) P a ,h{u) — (h 2 A — l)u, <VG$(P a ,„). 
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It is the self-adjoint operator associated to the semi-bounded quadratic form 
(i-iS) ^ 2 ||Viy||^ - ||u||m - ha\\4dM- 

The operator P a /, is linked with the semiclassical Dirichlet-to-Neumann 
operator as follows: if f is an eigenfunction of R sc |(A) with eigenvalue a, then 
the corresponding Helmholtz function u defined by (1.1) is in the domain 
(1.13) of P a h, and P 3i hU = 0 (where h = A -1 ). 

Then we have the following result, proved in Section 2. 

Proposition 1.1. Let h — A -1 . Assume 0 < 6-y < 0 2 < 27r. Then the 
number of eigenvalues of C(A) in the interval Je 1 ,e 2 '■= { e ‘ 6 ■ & G [ 61 , 62 )} is 
equal to 

(1.16) N(A; e 1 ,e 2 ) = N ( A; a 1: a 2 ) = Uf(a 2 ) - N^), 
where aj = — cot(#y/ 2 ) and 

(1.17) N h (a) := 7 \ \i is an eigenvalue of P a ^, /J < 0}. 

Having thus reduced the problem to a standard question about semiclas¬ 
sical spectral asymptotics, we obtain (after some calculations in Section 3) 
our main result. 

Theorem 1.2. (%) The following estimate for the quantity (1.11) holds: 

(1.18) N(A; a 1( a 2 ) = 0(A d_1 ); 


(ii) Further, if the set of periodic billiards on M has measure 0 then the 
following asymptotic holds as A —>• + 00 ; 

(1.19) N(A; a 1: a 2 ) = («(a 2 ) — K(a x )) yo\\dM)X d ~ 1 + o(A d_1 ), 

where K(a) is given explicitly by 


(1.20) «(a) 


Wd -1 
( 27 r ) d ~ 1 




_ rl 2,(d-l)/2_ 


a 2 + ry 


dg 


1 - + H(a)(l + a 2 r- 


l)/2 


Here H(-) is the Heaviside function, ojd is the volume of the unit ball in 
R d , and vol(M) and vol \ 8 M) are d-dimensional volume of M and (d — 1)- 
dimensional volume of dM respectively. 
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(in) In the case d — 3, we can evaluate this integral exactly and we find that 


( 1 . 21 ) 


«( a ) 


(-7 — - arccot(a)(l + a 2 ) + (1 + a 2 ) + -a) 
4n \ 4 7r 7 r / 


where arccot has range ( 0 , 7 r). This is simpler expressed in terms of 6. 
Defining k{6) — n(a) where a = —cot(0/2) = cot( 7 r — 9/2), we have (still 
under the zero-measure assumption on periodic billiards) 


( 1 . 22 ) 

(1.23) 


N(A; 9 lt 0 2 ) = (k(0 2 ) - k( 0!)) vol'(<9M)A 2 + o(A 2 ), 

~/a\_ 1 ( 1 1 / 9 - sin 9 

47 T V 4 + 27T ^sin 2 (0/2) 


K 


Remark 1.3. (i) It looks disheartening that the remainder is just “o” in 
comparison with the principal part and only under geometric condition of 
the global nature but it is the nature of the beast. Consider M a hemisphere; 
then for A 2 = n(n + d — 1) with n E Z + the operator R( A) has eigenvalue 0 
of multiplicity x n d ~ 1 and therefore we do not have even an asymptotic but 
just an estimate. 


(ii) We can explain this by pointing out that the problem we consider is 
truly d-dimensional, in contrast, for example, to the problem of distribution 
of negative eigenvalues of H — A with 

£(/-/) = {u 6 H 2 (X ) : d u u\ 9M + a/\^ 2 u\ dM = 0} 


with At, positive Laplacian on dM and a > 1 which is in fact (d — 1)- 
dimensional. 


(iii) On the other hand, calculations of Section 4 imply that for a / 0 
quantum periodicity will be broken (see [Ivr2], Section 8.3) and one can get 
remainder estimate o(A d_1 ) albeit with the oscillating principal part. 

(iv) Under certain assumptions about geodesic billiard flow remainder (see 
[Ivr2], Section 7.4) estimate (3.2) could be improved to 0(h 1 ~ d /\ log h\) 
or even 0(h 1 ~ d+s ) with a small exponent d > 0 and therefore remainder 
estimate (1.19) could be improved to 0(A d ^ 1 / log A) or even 0(A rf-1 ~' 5 ). 

(v) Despite different expressions for a < 0 and a > 0 one can observe that 
ft(a) is monotone increasing and continuous. 
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(vi) In particular, k( 0) = |(27r)” d a; c /_i and n(—oo) = — ^(2ir)~ d Ud-i. Since 
a — 0 corresponds to the Neumann boundary condition and a = — oo to 
Dirichlet (see Proposition 4.1), this agrees with (1.4). 

(vii) One can consider eigenvalues of operator p\~ 1 R(\) with p > 0 smooth 
on dM] then estimates (3.1), (1.18) and asymptotics (3.2), (1.19) hold in 
the frameworks of Statement (i) and (ii) of Theorem 1.2 respectively albeit 
with n{a) \jo\\3M) replaced by 



n(p(x')a) do 


where do is a natural measure on dM; however without this condition p > 0 
problem may be much more challenging; even self-adjointness is by no means 
guaranteed. 


(viii) Operators of the form P a ^ were considered by Frank and Geisinger 
[FG]. They showed that the trace of the negative part of P a h has a two-term 
expansion as h —y 0 regardless of dynamical assumptions 1 , and the second 
term in their expansion (the term of [FG, Theorem 1.1]) is closely related 
to hi(a) — see Remark 3.4. 

Remark 1.4. We can rephrase Theorem 1.2 in terms of a limiting measure 
on the unit circle. For each A > o, let p(h), h = A 1 , denote the atomic 
measure determined by the spectrum of C(A): 


(1.24) n(h) = {2nh)->- 1 Y, 

e' B J £specC(/?- 1 ) 


where we include each eigenvalue according to its multiplicity as usual. 
Then Theorem 1.2 can be expressed in the following way: the measures p(h) 
converge in the weak-* topology as h —y 0 to the measure 

(1.25) uJd-iMo\'(dM)-^-k{6)d6 on (0,27 t), that is on S 1 \ {1}. 

du 

In particular, this measure is absolutely continuous, and finite away from 
e' 9 = 1 with an infinite accumulation of mass as 0 f 2tt. In this form, we can 

1 The fact that Frank and Geisinger obtain a second term regardless of dynamical 
assumptions is simply due to the fact that they study Tr f(P a ,h) with f( A) = —XH(—X) 
( H is the Heaviside function), which is one order smoother than f(A) = H{— A). 
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compare our result with results on the semiclassical spectral asymptotics 
of scattering matrices. In [DGHH] and [GHZ], the scattering matrix S/,(E) 
at energy E for the Schrodinger operator h 2 A + \/(x) on was studied in 
the semiclassical limit h —> 0. Assuming that \/ is smooth and compactly 
supported, that E is a nontrapping energy level, and that the set of periodic 
trajectories of the classical scattering transformation on T*S d ~ 1 has measure 
zero, it was shown that the measure p(h) defined by (1.24) converged weak-* 
to a uniform measure on S 1 \ {1}, with an atom of infinite mass at the 
point 1. On the other hand, for polynomially decaying potentials, it was 
shown by Sobolev and Yafaev [SoYa] in the case of central potentials and 
by Gell-Redman and the first author more generally (work in progress) that 
there is a limiting measure which is nonuniform, and is qualitatively similar 
to the measure for C(/)" 1 ) above in that it is finite away from 1, with an 
infinite accumulation of mass at 1 from one side. 


2 Reduction to semiclassical spectral 
asymptotics 

In this section we prove Proposition 1.1. This result actually follows directly 
from the Birman-Schwinger principle. As some readers may not be familiar 
with this, we give the details. 

Proof of Proposition 1.1. We begin by recalling that the operator P a h is the 
self-adjoint operator associated to the quadratic form (1.15), that is, 

Q*M := h 2 ||Vu||^ - \\u\\ 2 m - ha\\u\\ 2 dM . 

We recall the min-max characterization of eigenvalues: the nth eigenvalue 
p n {a, h ) of P a h is equal to the infimum of 

sup Q a ,h(v) 

v£V,\\v\\=l 

over all subspaces V G H 1 (M) of dimension n. The monotonicity of Q a ^ 
in a, for fixed h, shows that the eigenvalues are monotone nonincreasing 
with a. In fact, they are strictly decreasing, which follows from the fact- 
that eigenfunctions of P a h cannot vanish at the boundary. Indeed, the 
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Figure 1: Diagram showing the variation of eigenvalues p(a, h) of P a /, as 
a function of a for fixed h. The eigenvalues are strictly decreasing in a. 
Consequently, the number of negative eigenvalues of P a2 ,h is equal to the 
number of negative eigenvalues of P ai ,h together with the number that cross 
the a- axis between a — a 1 and a — a 2 . 


eigenfunctions satisfy the boundary condition hd u u — — au , which shows 
that if u vanishes at the boundary, so does d u u, which is impossible. 

The eigenvalues /i„(a, h ) are thus continuous, strictly decreasing functions 
of a. Let ai < a 2 be real numbers. The Birman-Schwinger principle [Ivrl, 
Prop. 9.2.7] says that the number of negative eigenvalues of P 32 ,h is equal to 
the number of negative eigenvalues of P ai ,h, plus the number of eigenvalues 
/i„(a, h ) of P a ,h that change from nonnegative to negative as a varies from ai 
to a 2 . A diagram makes this clear: see Figure 1. 

The strict monotonicity of /i(a, h) in a shows that the number of eigen¬ 
values p n (a, h ) of P a h that change from nonnegative to negative as a varies 
from ai to a 2 is the same as the number of p(a, h ) (counted with multiplicity) 
equal to zero, for a £ [ai, a 2 ). Next, we observe that the space of eigenfunc¬ 
tions u„(a, h) of P a h with zero eigenvalue, i.e. p n (a, h) — 0 is in one-to-one 
correspondence with the space of eigenfunctions of C(A), A = h -1 , with 
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eigenvalue (a — i)(a + /) 1 , or equivalently e' e where a — — cot(0/2). Indeed, 
whenever u r is such an eigenfunction of P a then 

(2-1) f := iu)\ QM 

is an eigenfunction of C(A), with eigenvalue (a — /)(a + /) _1 . Conversely, 
if f is an eigenfunction of C(A) with eigenvalue (a — /)(a + /) _1 , then by 
definition there exists a Helmholtz function u such that u is related to f 
according to (2.1), and we have (hd u + a)u — 0 at 8M. This completes the 
proof. □ 

Remark 2.1. We can apply similar arguments for A -< 5 p - 1 /?(A) as p > 0 is 
a smooth function on dM and then plug corresponding parameters in the 
boundary conditions coming again to equality (1.16). 

We next digress to prove that the eigenvalues of C(A) are monotonic 
(that is, they move monotonically around the unit circle) in A. This plays 
no role in the remainder of our proof, but is (in the authors’ opinion) of 
independent interest. 

Proposition 2.2. The eigenvalues of C{ A) rotate clockwise around the unit 
circle as A increases. 

Remark 2.3. This implies that the eigenvalues of R sc |(A) are monotone 
decreasing in A. 

Proof. As discussed in the previous proof, C(A) has eigenvalue e 10 if and 
only if P a h has a zero eigenvalue, where a — a(9) — — cot(0/2). Thus, as a 
function of h — A -1 , 9(h) is defined implicitly by the condition 

p(a(9), h) = 0. 

Since a is a strictly increasing function of 9, and we have just seen that p is 
a strictly decreasing function of a, it suffices to show that when p = 0 , p is 
a strictly increasing function of h, hence a strictly decreasing function of A. 

We now compute the derivative of p with respect to h, at a value of a 
and h where p(a, h) — 0. We have 

-^p{a, h ) = ^((^ 2 A - 1 )u(h), u(h)) M 
= 2h(Au, u) M + ((h 2 A - 1 )u'(h), u(h)) M + ((h 2 A - l)u(h), u'(h)) M 
— 2h(Au, u) M + ((h 2 A - 1 )u, u) M ~ W, (h 2 A - 1 )u) M - 
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In the third line, we used the fact that (h 2 A — 1 )u — 0 when /i(ti) = 0. Note 
the second term is not zero, as u' is not in the domain of the operator due 
to the changing boundary condition, so we cannot move the operator to the 
right hand side of the inner product without incurring boundary terms. We 
use the Gauss-Green formula to express the last two terms as a boundary 
integral: 


p{h) = 2h(Au, u) M + h(hd u u, u) dM - h(u', hd v u) 9M 
= 2h(Au, u) M + h(hd u u\ u) dM + ha(u', u) dM . 

Differentiating the boundary condition we find that 

( hd u + a)u' — —hd v u at dM. 

Substituting that in we get 

p'(h) = 2h(Au, u) M - h(u, d u u) 9M - 

Applying Gauss-Green again, we get 

= h(Au, u) M + h\\Vu\\ 2 M 
= h 1 (\\u\\ 2 L 2( M ) + ||Wu|| 2 , 2(m) ) > 0. 


□ 


3 Semiclassical spectral asymptotics 

In this section, we prove Theorem 1.2. Essentially, we have arrived at a 
rather standard semiclassical spectral asymptotics problem and results are 
due to [Ivrl], Chapter 5 or [Ivr2], Chapter 7. 

Proposition 3.1. (i) Let N h (a) be as in (1.17). The following asymptotic 
holds as h —>■ +0: 


(3.1) 


N-(a) = (2nh)- d u d m\(M) + 0(11 1 -") 
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(ii) Further, if the set of periodic billiards on M has measure 0 then as 

h -> +0: 

(3.2) N^(a) — (2nh)~ d Ud vol(M) + h^ d n(a)\/o\'(dM) + o(h 1 ^ d ) 
with n(a) given by (1.20). 

Proof. One can check easily that the operator P a j, is microhyperbolic at 
energy level 0 at each point (x, £) G T* M in the direction £; further, the 
boundary value problem is microhyperbolic at each point (x'; £') G T*dM 
at energy level 0 in the multidirection (£', ) with £i = roots of 

— 0 ; finally, the boundary value problem is elliptic at each point 
of the elliptic zone (c T*dM ) if a < 0 , and either elliptic or microhyperbolic 
in the direction £' at each point of the elliptic zone (c T*dM ) if a > 0 - 
see definitions in Chapters 2, 3 of [Ivr2], Then statements (1.18), (1.19) 
follow from Theorems 7.3.11 and 7.4.1 of [Ivr2], 

We now assume that the set of periodic billiards on M has measure zero, 
and compute the second term in the spectral asymptotic explicitly. Similar 
calculations appear in [FG]. 

To do this, one can use method of freezing coefficients (see f.e. [Ivr2], 
7.2) which results in 

(3.3) h 1 ~ d K,(a) = f (e(0, xp 0, xp 1) — (2nh)~ d ujd) dxi 

J E+ 

where e(x', xp y', yp r) is the Schwartz kernel of the spectral projector E{t) 
of the operator H a — h 2 A in half-space M d_1 x M + 9 (x'.xx) with domain 
D(H a ) = {u G H 2 : (hd Xl + a)u\ Xl=0 = 0}. 

We obtain this spectral projector by integrating the spectral measure. 
This in turn is obtained via Stone’s formula 

(3.4) dE L (o) = ((/_ - (a + /0 ))" 1 — (L - (a — /0) -1 ) da. 

Consider the resolvent for H a , (H a — a) 1 , for a G C \ M. Using the 
Fourier transform in the x' variables, we can write the Schwartz kernel of 
this resolvent in the form 

(2>r hf~ d J e'<''T a + |{'| 2 - ctJ'V, n) d?. 


(3.5) 
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Here T a is the one-dimensional operator T a = —h 2 d 2 + |£'| 2 on /. 2 (M + ) under 
the boundary condition (hd + a)u| Xl= o = 0. The spectral projector Eh 3 ( 1) is 
therefore given by 

(3.6) (2nhy- d I I e l{x '- y ' H 'dE Ta {(r - |£'| 2 )( x i, yi) df do. 

J — OO J 

Thus, we need to find the spectral measure for T a . Write o — |£'| 2 = p 2 , 
where we take p to be in the first quadrant of C for Im o > 0, and in the 
fourth quadrant for Imcr < 0. 


Lemma 3.2. Suppose that \mp > 0 and Rep > 0. Then the resolvent kernel 
( T a — r/ 2 )” 1 takes the form 


(3.7) 

If Im p 
form 


C T a -p 2 ){x,y ) 


_L_ f e ir/(x-y)/h iy-a c iy(x+y)/h\ 

2/777 V /77+a j 1 

J— ( e iv{y-x)/h 1 iy-a c m(x+y)/h\ 
2/777 \ /77+a J 1 


< 0 and Rep > 0, then the resolvent kernel (T a — 


x > y 
x < y. 

p 2 ) _1 takes the 


(3.8) ( T a — p 2 )(x, y) 


_ L( e iy(y~x)/h 1 i’l+a e -iy(x+y)/h\ 

2/777 \ irj—a J ’ 

e iy( x ~y)/ h 4 . in ±j e -hh+r)/^ 

2 hr) V irj—a I ’ 


x > y 
x < y. 


Proof. In the regions x < y and x > y. the resolvent kernel must be a linear 
combination of e" ,X///l and e~ IT,x l h . Moreover, for lm /7 > 0, we can only have 
the e +lvx ^ h term, as x —» 00 , as the other would be exponentially increasing. 
So we can write the kernel in the form 


(3.9) 


eye 


y-iijx/ h 


x > y 


c 2 e +irix/h + c 3 e~ irix/h , 


x < y. 


We apply the boundary condition, and the two connection conditions at 
x = y, namely continuity, and a jump in the derivative of — 1/h, in order to 
obtain the delta function h(x — y) after applying T a . These three conditions 
determine the c ( - uniquely, and we find that, in the case Im p > 0 , 


(3.10)! 

c 1= ' ( e ~iyy/h + n l a e +iyy/h\ 

2 hp v ip + a ’ 

(3.10)2,3 

- _ ' ' r l ~ a „+iyy/h 

^2 a/ ■ | ^ » ^3 r 

Zhr] it] + a ^ 


-iyy/h 


yielding (3.7). A similar calculation yields (3.8). 


□ 
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We now apply (3.4) to find the Schwartz kernel of the spectral measure 
for T a . 

Lemma 3.3. The spectral measure dEjX T ) is given by the following. 

(i) For t > 0, r — rj 1 

(3.11) dE T ,( t) 

. f c iy{x-y)/h | c iri(y-x)/h | lr 1 ~ 3 c iy(x+y)/h | 'V + a c -ir)(x+y)/h\ 2tjdq 

V in + a in — a J 


An h r] 


(ii) For t < 0, the spectral measure dE(r) vanishes for a < 0 , while for 

a > 0 


(3.12) 


cIEtXt) = + a 2 )dr. 


Proof. This follows directly from Lemma 3.2 and Stone’s formula, (3.4). 
The extra term for a > 0 arises from the pole in the denominator, iq + a 
for Im q > 0 and iq — a for Im 77 < 0 in the expressions (3.7), (3.8), which 
only occurs for a > 0. For r negative, we need to set 77 = i\f—r + 0 in 
(3.7) and 77 = —i\f—r + 0 in (3.8), and subtract. Then everything cancels 
except at the pole, where we obtain a delta function —2ni8{y/~ T — a), which 
arises from (yf—r + /0 + a) -1 — (y/—T ~ /0 + a) -1 . This term arises from 
the negative eigenvalue —a 2 which occurs for a > 0, corresponding to the 
eigenfunction \j2a/h e~ ax / h . □ 

Plugging this into (3.6), and making use of the fact that dodf — 2qdqdf', 
we find that the Schwartz kernel of Eh . (1) is given by 

(3.13) (2nh)~ d J I H{ 1 - |£'| 2 - q 2 )e i{x '- y ' H ' 

x ( e ir i(xi-yi)/ h _|_ e iv(yi-xi)/h _|_ XL ~ a c m(x\+yi)/h | ^_Zt_^ e -''?( x i+yi)//A 

V iq + a iq — a / 

for a < 0 while for a > 0, it is given by the sum of (3.13) and 


/ 2 a 


1 — r -axy!h r -ay\/h 

h 


(3.14) ( 2nhy “ — e 


Ax'-y'K' 


5(<x - l£l 2 + a 2 )h£' da. 
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We are actually interested in the value on the diagonal. Setting x = y, and 
performing the trivial £' integral, we find that the Schwartz kernel of the 
spectral projector E Ha ( l)(x, x) on the diagonal is given by 


(3.15) TWhW / (! - r/ 2 ) (d ~ 1)/2 2 + ——-e 2/,?Xl//? + df) 


Ud- 


(27 rh) c 


H(a) 


n 7 + a 

(d - l) co d -i a n - 
(2nh) 


rr] — a 


d— 


- 2 ^/" / (n + a 2 )( d - 3 )/ 2 dn. 


Since 




_! / 2(1 — 77 2 )( d-1 )/ 2 dr ] — 


OJd, 


we see by comparing with (3.3) that this term disappears in the expression 
for /t(a) and we have, after performing the x\ integral as in (3.3) 

. _ , fiii * r 1 

(3.16) h 1 d «(a) = 


(2nh) c 


( l- ?/ 2 )(a-i)/2 


x 


h ~ af %+ior' 


i rj + a V 2 


ir) + a(iK_ /Q ) i\ | drj 


irj — a V 2 


H(a) 


(d - l)cu rf _i f 1 


2(2irh) d 1 J_ 

Simplifying a bit, and performing the integral, we have 


(cr + a 2 )( d_3 )/2 da 


(3 ' 17) K(a) = /V -» 2 ) l ‘ , - ,)/2 ^ T ^(» + - 0)- 1 */ 

+ H(a) (Sph (1 + a2)( ‘'" 1,/2 ' 

We further simplify this expression by expanding ( irj — a) 2 — a 2 — 2/a// — r/ 2 , 
and noting that the contribution of the — rj 2 term is zero, as this gives an 
odd integrand in the // integral. A similar statement can be made for the a 2 
term, except that there is a contribution from the pole in this case. This 
leads to the expression 


(3.18) «(a) 


Ud -1 
( 27 r ) d ^ 1 



drj 


1 

4 


H(a)(l + a 2 )< d - 


l)/2 
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Although not immediately apparent, this formula is continuous at a = 0. In 
fact, the function a(a 2 + r/ 2 ) -1 has a distributional limit (sgn a) TrS(rj) as a 
tends to zero from above or below. The change of sign as a crosses 0 means 
that the integral in (3.18) has a jump of —1 as a crosses zero from negative 
to positive. That exactly compensates the jump in the final term. 

In odd dimensions, we can compute this integral exactly. In particular, 
in dimension d — 3, we find that 


(3.19) 


K ( a ) 


0J2 

(2vr) 2 


i + ^ + (i + a 2 )(i- 


arccot a 

7r 



□ 

Proof of Theorem 1.2. This follows immediately from Proposition 3.1 and 
Proposition 1.1. □ 

Remark 3.4. The second term of the expansion in [FG, Theorem 1.1] is 
obtained by computing 

(3.20) (2vr hY~ d [ (1 -a) [ dE Ta {a - \e\ 2 ){x ltyi ) d? da 

J — OO J 

instead of (3.6). 


4 Relation to Dirichlet boundary condition 

In this section we observe that the limit a —> —00 corresponds to the 
Dirichlet boundary condition. More precisely, we have 

Proposition 4.1. Let N^(— 00) denote the limit 

N^(— 00) := lim Nj((a), 

a—>— 00 

where (a) is given by (1.17). Then we have 

(4.1) N ft (— 00) = #{Aj < /P 1 \ is a Dirichlet eigenvalue of A}. 

Remark 4.2. Because the quadratic form (1.15) is monotone in a, the counting 
function N^(a) is monotone in a. Hence the limit above exists. 
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Proof. We use the rnin-rnax characterisation of eigenvalues. Let No(A) 
denote the number of Dirichlet eigenvalues (counted with multiplicity) less 
than or equal to A = h _1 . This is equal to the maximal dimension of a 
subspace of /4 q(M) on which the quadratic form Q D , given by 

(4.2) Q D (u,u) = h 2 \\Vu\\ 2 2 -\\u\\ 2 

is negative semidefinite. On the other hand. Nf (a) is equal to the maximal 
dimension of a subspace of H 1 (M) on which the quadratic form Q a given by 
(1.15) is (strictly) negative definite. 

We first show that No(fj 1 ) < (—oo). Let V be the vector space 

spanned by Dirichlet eigenfunctions with eigenvalue < A 2 . Clearly, the 
quadratic form Q a is negative semz definite on \/, and if A 2 is not a Dirichlet 
eigenvalue, then it is negative definite, proving the assertion. In the case that 
A 2 is a Dirichlet eigenvalue, we perturb V to V e , a vector space of H 1 (M) of 
the same dimension as \/, so that, for for e sufficiently small depending on 
a, Q a is negative definite on V € . For simplicity we only do this in the case 
that the A 2 -eigenspace is one dimensional, leaving the general case to the 
reader. To do this, we choose an orthonormal basis of V (with respect to 
the L 2 inner product) of Dirichlet eigenfunctions v 1 ,... ,v k with eigenvalues 
A 2 ... A 2 , where \ k — A. Then we perturb only v k , leaving the others fixed. 
We choose s G /-/q(/W) ± , the orthogonal complement of Hq(M) in H 1 (M) 
(with respect to the inner product in so that 

(4.3) Q a (y-, s) = 0, i < k and Q a (v k , s) > 0. 

We check that this is possible. Notice that s E Hq(M) 1 implies that 
(A + l)s = 0 in M. Then as v, has zero boundary data, we have 

(4.4) (A- + l)(v7, s)m = (A Vi,s) M - (vi, As) M = {d v v il s) dM . 

We choose s so that (d u \ /■,, s)g M vanishes for / < k and is positive for / = k. 
This is possible: in fact, due to the unique solvability of the boundary value 
problem 

(4.5) (A + l)s = 0, s| 9/w =fe H 1/2 (M), 

for s E we see that s can have any boundary value in H 1 ^ 2 (dM) 

which is dense in L 2 (dM). Then using (4.4) we see that (d u v k ,s)gM > 0 
implies that ( v kl s)M > 0. 
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We now define V e to be the span of v \,..., v k -\ and v k + es. Then we 
have 

(4.6) Q a {v it v k + es) = 0, i < k 
and 

(4.7) Q a {v k + es, v k + es) = Q a (v k , v k ) + 2 eQ a (v k ,s) + e 2 Q a (s,s) 

= 2eQ a (v il Sj) + e 2 Q a (s i ,Si) 

— —2e(h 2 + l)(v;, s;)m + 0(e 2 a 2 ) 

which is strictly negative for ea 2 small enough. It follows that Q a is negative 
definite on V k when ea 2 is small enough. A similar construction can be made 
when A 2 has multiplicity greater than 1. 

We next show that No(fi 1 ) > (—oo). We argue by contradiction: 
if not, then for any a, there is a vector space W of dimension > k + 1 
on which Q a is negative definite. Then there is a nonzero vector w G 1/1/ 
orthogonal (in the H 1 (M) inner product) to V. We can write w — w' + s 
where w' G Hq(M) and s G Hq(M) 1 . Then w' is a linear combination of 
Dirichlet eigenfunctions with eigenvalue > X' > X, where X' is the smallest 
eigenvalue larger than A. We then have 

(4.8) 0 > Q a (w’ T s, w’ T s) — Q a (w', w') + 2 Q a (w', s) + Q a (s, s) 

> (A' - X)\\w'\\l - 2 (h 2 + 1 )||m/|| 2 ||s|| L 2 (m) - ha||s|| 2 2(aM) . 

However, some standard potential theory shows that ||s||^ 2 ( M ) is bounded 
by a constant times ||s|| t 2 ( 9/W )- T° see this, extend M to a larger manifold 
M of the same dimension, and let G(x,y) be the Schwartz kernel of the 
inverse of (A ^ + 1) _1 on L 2 (M), with Dirichlet boundary conditions at dM. 
We can write s as f. )M d Uy G(x, y)h(y) dy where (1/2 + D)h — s\sm and D 
is the double layer operator on dM determined by G. Standard arguments 
show that (1/2 + D) has a bounded inverse on L 2 (dM) and d Uy G(x,y) is 
a bounded integral operator from L 2 (dM) to L 2 (M). So we can write, for 
a < 0, 

0 > Q a (w' + s, w’ + s) 

> (A' - A)|| la/'Hl - 2 C(h 2 + 1)|| w'|| 2 1|^||/. 2 (6»/w) + ^l a lll s llr 2 (£>M) 2 
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and the RHS is clearly positive for |a| large enough, giving us the desired 
contradiction. 

□ 
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